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ABSTRACT 



Knowledge of Jupiter's deep interior would provide unique constraints on the formation of the Solar System. 
Measurement of its core mass and global composition would shed light on whether the planet formed by accretion 
or by direct gravitational collapse. At present, the inner structure of Jupiter is poorly constrained and seismology, 
which consists of identifying acoustic eigenmodes, offers a way to directly measure its deep sound speed profile, and 
thus its physical properties. Seismology of Jupiter has been considered since the mid 1970s, but hitherto the vari- 
ous attempts to detect global modes led, at best, to ambiguous results. We report the detection of global modes of 
Jupiter, based on radial velocity measurements performed with the SYMPA Fourier spectro-imager. The global seismic 
parameters that we measure include the frequency of maximum amplitude 1213 =b 50 /iHz, the mean large frequency 
spacing between radial harmonics 155.3 ± 2.2 /iHz and the mode maximum amplitude 49^ cm s_1 ' an values that are 
consistent with current models of Jupiter. This result opens the way to the investigation of the inner structure of the 
Solar System's giant planets based on seismology techniques. 
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1. Introduction 

Jovian seismology has long been considered as both a po- 
tentially powerful to ol for probing the interior of Jupiter 



([Guillot et al.l 120041) and a natural extension of helio- 
seismology ([Goldreich fc Keelevl Il977| ) , because the com- 
mon fluid nature of Jupiter and the Sun is expected 
to lead to similar oscillations and the possibility of us- 
ing similar observational techniques. Theo retical works 
(jVorontsov et aLlll976t iBercovici fc Schubertlfl987l ) predict 

■ that Jovian global oscillations should have a frequency 
range of [800, 3500] /iHz with 10 to 100 cm s _1 amplitude, 
values that are comparable to those of the Sun. Several at- 

■ tempts to observe Jovian global modes have been carried 
out s ince the mid-1980s , with thermal infrared photom- 
etry ([Deming et al.lll989 | ). and Doppler spec trometry us- 

' ing both magneto-optic (Schmider et al. 1991) and Fourier 
transform spectrometers ([Mosser et al.|[l993ll2000f), respec- 
tively. Moreover, in frared observations ([Walter et al.l H 996; 
M osser et all 119961) were performed to look for pressure 
waves excited by the impact of the comet Shoemaker-Levy 
9. Infrared observations are affected by atmospheric in- 
homogeneities and have been unsuccessful so far. In con- 
trast, all of the Doppler measurements exhibit an excess 
power of about 1 m 2 s -2 between 800 and 2000 /iHz, 
which cannot be explained by instrumental systematics nor 
spurious atmospheric signals. A tentative comb-like struc- 
ture with a 139 ± 3 /iHz mean spacing was also identi- 



fied, but the combination of a low signal-to-noise ratio 
(hereafter SNR), this frequency corresponding to the least 
common multiple of terrestrial and Jovian rotation fre- 
quencies, and th e spacing being incompatible with theo- 
retical estimates (iProvost et al.lfl993b iGudkova et al.| [l995; 
iGudkova fc Zharkovlll999f ) implied that it was probably an 
artifact. The need for specific instrumentation able to com- 
bine high spectral and spatial resolution thus emerged from 
the 1980 and 1990s experiments. The major difficulty in 
performing seismic observations of Jupiter is related to its 
rapid rotation (500 m s _1 at the equator on 1 arcsec), which 
diminishes the instrument's velocity sensitivity and makes 
it extrem ely sensitive to pointi ng errors. The SYMPA in- 
strument (ISchmider et al.l l2007, hereafter S07) is a Fourier 
tachometer designed to circumvent this difficulty, whose 
principle is based on the spectro-imaging of the full plan- 
etary disk in a non-scanning mod e, and is related to the 
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helioseismic instr uments GONG ( Beckers fc Brownl Il979f ) 
and MDI/SOHO ([Scherrer et aDll995h . It produces radial 
velocity maps of Jupiter's upper troposphere by measur- 
ing the Doppler shift of solar Mg lines at 517 nm that are 
reflected by Jupiter's clouds (S07). 

Two instruments were built and simultaneously used 
during two observation runs in 2004 and 2005 at San Pedro 
Martir and Teide observatories. The best observation se- 
quence with SYMPA was acquired as part of the 10-day 
run performed at the Teide Observatory between April 2 
and 12 2005. Jupiter was then in opposition with positive 
declination. Such favorable conditions for observations did 
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not occur again in the northern hemisphere until January 
2011. The data were found to be of much higher quality 
than all previous observations: at the end of the run, with 
a 21.5% duty cycle, the mean noise level reached 12 cm 
s _1 per sample, which was five times lowe r than previously 
achieved for Jupiter (Mos ser et alj|2000l ). At the time, a 
prelimi nary analysis did no t yield the detection of Jovian 
modes ([Gaulme et alJl2008l herafter G08). 

Here we introduce the detailed analysis of the 2005 data 
obtained at Teide Observatory, where we focus on the time 
series with the highest SNR. A short review of the data 
processing explains the procedure that we followed to op- 
timize our chance of finding out the global modes, instead 
of performing a global analysis as in G08. This leads to the 
detection of an excess power modulated by a comb struc- 
ture, whose properties are compatible with Jovian modes 
(Sect. 2). We then explore all of the possible origins of the 
detected signal to conclude of the detection of Jupiter's 
global modes (Sect 3). Last, we conclude on the implica- 
tion that this result generates about future studies of giant 
planets (Sect. 4). 



2. Detailed analysis of the SYMPA data 

2.1. Data processing short review 

We review the essential steps of the whole data processing, 
from a single acquisition to the production of power spectra 
associated with a set of spherical harmonics. We refer the 
reader to S07 and G08 for a detailed description. 

SYMPA is a Fourier tachometer based on a Mach- 
Zehnder (MZ) interferometer working at fixed optical path 
difference. It produces four interferograms of the solar spec- 
trum reflected by Jupiter, in a 5-nm spectral range cen- 
tered around the Mg triplet at 517 nm. Each interfero- 
gram modulates the Jovian photometric figure. The dif- 
ferences between the two pairs of images allow us to can- 
cel the photometric intensity, to produce two fringe pat- 
terns that are phase-shifted by about 90 degrees. This step 
involves the application of an algorithm that straightens 
up each image from geometrical distortions, then overlaps 
the two pairs of images. However, the accuracy of this 
process was insufficient to fully cancel the Jovian image 
in both interferograms. In particular, Jupiter's zones and 
belts were still slightly visible. Thus, we applied a Wiener 
filter in the 2D Fourier transform of both interferograms 
in order to substantially reduce the photometric residu- 
als. Moreover, towards the limb, the SNR is low because 
of Jupiter's limb darkening, while geometrical distortions 
are maximum. Therefore, the velocity maps are limited to 
the central 75% of the Jovian radius, applying an apodisa- 
tion at the edge. 

The linear combination of the two interferograms leads, 
after substraction of the instrumental phase, to a radial 
velocity map. In this map, the Jovian rotation is the ma- 
jor effect. The first step consists of subtracting from the 
velocity map the Jovian rotation plus the relative motion 
of Jupiter with respect to the observer. This requires us 
to precisely determine the inclination of Jupiter's rotation 
axis on the image, which is achieved within a fraction of a 
degree thanks to the belts and zones on Jupiter. The ro- 
tation law is then applied to the data after centering the 
image. The precision of the centering is one of the limiting 
factors of the instrument's performance. The velocity maps 
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Fig. 1. Representation of the projected spherical harmon- 
ics, as a function of degree I and azimuthal order m, until 
degree i = 5. The masks used to filter the velocity maps 
are limited to 75% of the Jovian diameter. 

Table 1. Coupling between the masks, with value larger 
than 20%. 



£ = 
m = 


£ = 1 £ = 2 £ = 2 £ = 3 £ = 3 
m = — 1 ra = m — 2 m — — 3 m = — 1 
98.0% -66.7% -79.5% -40.1% -33.8% 


£ = l 
m = 


£ = 2 £ = 3 £ = 3 £ = 4 £ = 4 
m — — 1 rn = m — 2 m — —3 m — —1 
97.0% -72.1% -77.0% -37.8% -40.7% 


£= 1 
m—\ 


£ = 2 £ = 3 £ = 3 £ = 4 £ = 4 
m — —2 m—\ m — 3 m— —4 rn — —2 
97.0% -45.2% -79.8% 45.4% 23.7% 


£ = 2 
m — 1 


£ = 3 £ = 4 £ = 4 £ = b £ = 5 
m — —2 rn — 1 m — 3 m — —4 m — —2 
-96.4% -56.6% -78.2% 43.7% 31.8% 



were then decomposed into spherical harmonics: each map 
was multiplied by a set of masks, made of the projections 
of spherical harmonics in two dimensions. All pixels of each 
filtered map were then weighted with spherical harmonics 
to create a set of time series. The power spectrum of each 
time series was computed with a discrete Fourier transform. 
Spherical harmonics were considered until degree £ = 25, 
as it matches the maximum spatial resolution of SYMPA. 

2.2. Data with best SNR 

The rapid global analysis of the 2005 dataset (G08) did not 
yield a detection of global oscillations of Jupiter because of 
three independent factors. 

Firstly, the calculation of the power spectra discovered 
a spurious signal, at 5920 ± 200 /iHz, which is strongly cor- 
related with the position of Jupiter on the detector. This 
signal is due to the imperfect disentangling of the photo- 
metric and spectrometric contents. In particular, its ampli- 
tude is mostly a function of the fluctuations in the position 
of Jupiter along the x-axis. This dependence on the posi- 
tion of Jupiter ensures that the time series associated with 
symmetrical masks is more sensitive to the systematic er- 
rors. Moreover, an additional spurious signal was found, 
with a central frequency proportional to the product of 
the azimuthal order m times the Jovian rotation frequency. 
From order m = 5, this spurious signal contaminates all 
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Fig. 2. Power spectrum of the time series associated with the spherical harmonics of degree 1=1. The thin grey and 
bold black lines represent the power spectrum smoothed with a moving average, weighted by a triangular window, with 
two different scales. The grey line uses a 3.1-/iHz wide window, and the black line a linearly increasing width window, 
from 14 to 129 Hz across the [500, 4500] /iHz frequency range. The red dashed line indicates the background level fit 
and the full red line the background plus excess power fit. The upper plot represents the power spectrum of the window 
function, normalized such that its maximum value equals 1. 



power spectra at frequencies higher than 800 /iHz. This ar- 
tifact impedes the search for Jovian global oscillations in 
the power spectra associated with spherical harmonics of 
azimuthal order m < 5. 

The second factor is related to the blurring of the Jovian 
disk that is generated by both atmospheric seeing and ran- 
dom guiding fluctuations. The mean standard deviation of 
the position of Jupiter on the detector was measured to be 
3.3 pixel rms, for the whole observation run, i.e, 2.1 arc- 
sec. Since the apparent diameter of Jupiter was 44 arcsec 
in April 2005, only 20 elements of resolution can be consid- 
ered to be independent along its diameter. Moreover, the 
seeing varied during the observations and could sometimes 
reach a value of 2 arcsec, which contributed significantly to 
the blurring. Because we limited the analysis to 75 % of the 
Jovian radius, this blurring prevents us from searching for 
modes of degree higher than ~ 8. It also signifies that the 
mode sensitivity decreases as the degree increases. 

Thirdly, the projection of spherical harmonics in two 
dimensions generates a non-ort honor mal base. The elim- 
ination of the external ring of 25% of Jupiter's diameter, 
which contributes more to the noise than to the signal, dra- 
matically increases the degeneracy of the projected base. In 
Table [TJ we present the coupling of these four masks with 
their corresponding correlated masks. It turns out that if we 
consider the spherical harmonics until degree 1=5, only four 
masks remain independent: (l,m = 0,0), (l,m = 1,0), (l,m 
= 1,±1), and (l,m = 2,±1). Instead of globally analyzing 
all of the available data, as performed in G08, the present 
analysis focuses on the time series with the highest SNR. 
The probability of finding modes in our data is maximized 
by dealing with: i) antisymmetric masks, ii) azimuthal or- 
ders m < 5, hi) degrees as low as possible, and iv) inde- 
pendent masks to avoid redundancy. Therefore, we here- 
after consider the two time series that are filtered with the 
masks associated with spherical harmonics of degree and 
azimuthal order: l,m =1,0 a nd l,m = 1,± 1. Unlike a classi- 
cal helioseismology analysis ([Schoulll992l ), we do not apply 
the Legendre polynomial function to each m value, but ap- 
ply only the antisymmetric masks (Fig. [J). Hence, the time 



series contains both prograde and retrograde modes when 
the azimuthal order is not null. This procedure allows us to 
eliminate most of the spurious signal that would globally 
affect the velocity map, as guiding error and MZ tempera- 
ture variations. 



2.3. A clear excess power modulated by a comb-like structure 

We sum the information contained in these two indepen- 
dent time series by computing the power spectrum of the 
mean time series. The data set constitutes 27946 measure- 
ments acquired between April 2 at 23:48:14 (UT) and April 
12 at 06:30:17 (UT), with a 6-s time sampling and a 21.5% 
duty cycle. The power spectrum is computed with the dis- 
crete Fourier transform with an oversampling of a factor 
of two and a frequency resolution of 1.248 /iHz (Fig. [2j). 
We then analyze the power spectrum with methods sim- 
ilar to those developed for helioseismology (|Harvevlll985f ) 
and successfully ap plied to astero seismic data from C oRoT 
([Michel et al.l l2QQ8h and Kepler (jChaplin et al.ll2010D . As 
in previous Doppler measurements, the time series exhibits 
excess power between 800 and 2000 /iHz and a secondary 
excess between 2400 and 3400 /iHz. This time it is modu- 
lated by a clear comb-like structure of high SNR. 

We fit the mean background with a semi-Lorentzian 
function and each power excess with a Gaussian function, 
using the maximum likelihood estimator by taking into 
account the nature of the statistics, a x 2 w ^h two de- 
grees of freedom. From this fitting, the frequencies of max- 
imum amplitude of both excess powers are 1213 ± 50 and 
2781 ± 49 /iHz, with amplitudes of 9.8 ± 3.0 cm s" 1 and 
5.0±2.3 cm s _1 , respectively. The white noise has an ampli- 
tude of 7. 7 ±1.0 cm s _1 . To estimate the splitting frequency 
of the comb-like structure, we calculate the Fourier trans- 
form of the power spectrum by restricting the frequency 
domain to the main excess power [800, 2100] /iHz. The re- 
sulting power spectrum as a function of time (Fig. [3j) is 
equivalent to an autocorrelation of the time series. It is 
characterized by a main peak at 1.780 h, with second and 
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Table 2. Frequencies and amplitudes of the peaks detected 
by the HO statistical test with 1% significance level. 



V 


Velocity 


Error 


V 


Velocity 


Error 


fjJiz 


cm s _1 


cm s _1 




cm s _1 


cm s _1 


792 


44.0 


-6.2/+3.9 


1478 


46.4 


-6.5/+4.1 


854 


46.7 


-6.6/+4.2 


1533 


37.3 


-5.3/+3.3 


915 


34.1 


-4.8/+3.0 


1615 


40.9 


-5.8/+3.7 


970 


48.7 


-6.9/+4.4 


1753 


33.0 


-4.6/+2.9 


1011 


51.4 


-7.2/+4.6 


1939 


32.0 


-4.4/+2.8 


1066 


45.7 


-6.4/+4.1 


2110 


30.1 


-4.2/+2.7 


1094 


42.4 


-6.0/+3.8 


2535 


30.3 


-4.3/+2.7 


1162 


54.1 


-7.6/+4.8 


2714 


30.6 


-4.3/+2.7 


1245 


53.8 


-7.6/+4.8 


2837 


36.2 


-5.1/+3.2 


1341 


51.5 


-7.3/+4.6 


2947 


41.1 


-5.8/+3.7 


1410 


40.7 


-5.7/+3.6 


3071 


30.7 


-4.3/+2.7 



third harmonics at 3.584 h and 5.373 h and a half harmonic 
at 0.845 h. A linear fitting of the four values yields a corre- 
lation time of 1.789 ±0.026 h, corresponding to a frequency 
spacing of A^o = 155.3 ±2.2 /iHz. The fact that we observe 
a half harmonic in time rather than only integer harmon- 
ics can be due either to a combination of the window ob- 
serving function with the modulation of nearby oscillation 
modes or to a stochastic excitation of the modes. We also 
detect the presence of a peak at half the Jovian rotation pe- 
riod as would be expected from the rotational splitting of 
1=1 azimuthal modes. The echelle diagram associated with 
this mean spacing frequency A^o (Fig.[5j) clearly shows two 
maxima in the ~ 1000- 1600 /iHz range that are separated 
by half of the 155 /iHz folding frequency. Furthermore, two 
secondary maxima of smaller amplitude mirror them in the 
2500-3000 /iHz domain. 

To assess the significance of the observed signal, we then 
apply a null hypothesis (HO) statistical test to the peaks 
identified in the power spectrum. We retain only peaks with 
a significance level of 1% and remove all neighboring aliases 
caused by the Earth's rotational frequency (±11.6 /iHz). 
This results in the identification of 22 distinct local max- 
ima between 800 and 3100 /iHz. Both amplitude and fre- 
quencies were measured from the power spectrum rebinned 
over 11 points, to smooth out the daily aliases. To esti- 
mate the true amplitude of the detected peaks, we intro- 
duced 40 pure sine curves with amplitudes corresponding 
to a Doppler shift of 1 m s _1 in the time series, at fre- 
quencies away from the excess power. Because of the duty 
cycle and the rebinning, the sine curves appear as peaks 
of mean amplitude of [12. 2^;^ cm s_1 ] 2 m the power spec- 
trum. This scaling factor was used to infer the peak heights. 
If we assume these peaks are of Jovian origin, we then have 
to divide each peak height by a factor two, to take into 
account the doubling of the Doppler effect that affects the 
shift of reflected lines. By supposing these peaks correspond 
to pure sinusoidal signals, the maximum velocity for these 
is 49l^o cm s_1 - Moreover, the frequency measurement of 
40 artificial sine curves in the data allowed us to estimate 
the error in the frequency to be -5/±9 /iHz. 

3. Evidence of Jovian p-modes 

Our measurements agree with theoretical expectations in 
terms of the fre quency range, the amplitude , and the mean 
large spacing ([Bercovici fc Schubertl Il987t iProvost et al.l 
Il993[ ) and would correspond to the signature of modes of 
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Fig. 3. Power spectrum of the [750-2050 /iHz] power spec- 
trum of the time series. It is equivalent to the square of the 
envelope of the autocorrelation of the time series |^4C| 2 , 
and is computed as the Fourier transform of the power 
spectrum filtered with a tapered cosine window centered 
on 1400 /iHz and with 650 /iHz half-window. It is normal- 
ized such that |AC(0)| 2 = 1, then filtered for correlation 
times lower than 0.3 h. On one hand, 4 peaks, labeled 
(1/2,1,2,3) A^q -1 , are associated with the mean spacing 
of the modes: Avo = 155.3 ± 2.2 /iHz. The correspond- 
ing period, 1.789 ± 0.026 h, matches twice the sound speed 
transit time between the centre and the surface of Jupiter. 
On the other hand, the second highest peak at 5.030 h 
almost matches twice the Jovian rotation period ^j up , a 
spacing that is compatible with the mode splitting between 
positive and negative azimuthal orders (e.g. l,m=l,±l and 
l,m=l,-l). 

consecutive degrees (e.g. 1=1 a nd 1=2). Oscillating stars 
also display simila r signatures ([Appourchaux et all 120081 : 
iGarcfa et al.l 120091 ). However, other effects could a priori 
cause the measured signal: temperature fluctuations in the 
interferometer, telescope pointing errors, albedo features 
on Jupiter, and solar p-modes reflected by Jupiter's atmo- 
sphere. 

3.1. The contribution from the instrumental systematics 

Two sources of noise were identified in the data (G08): i) 
Slow temperature fluctuations of the interferometer cause 
a spurious signal at frequencies lower than 800 /iHz. ii) 
Periodic variations in the position of Jupiter on the de- 
tector generate noise whose characteristic frequency is 
5920 ±200 /iHz, i.e., a period of 2.8 min. These two sources 
of spurious signals could a priori contain harmonics possi- 
bly explaining the measured excess power. To assess their 
contribution to the signal in the frequency range where we 
detect Jupiter's oscillation modes, in Fig. [6] we compare the 
autocorrelation of the power spectrum in three frequency 
ranges: the first [10, 500] /iHz is dominated by tempera- 
ture effects, the second [800, 3800] /iHz is where the power 
excesses are detected, and the third is [4450, 7350] /iHz. 
The autocorrelation of the first frequency range does not 
have any particular characteristics, as would be expected 
for a thermal fluctuation spectrum. Its time constant is 
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Fig. 4. Echelle diagram of the power spectrum for a fold- 
ing frequency Auq = 155 /iHz. The power spectrum was 
divided by the fitted background function and smoothed 
over five bins, i.e., 3.1/iHz. It was then divided into con- 
secutive intervals of 155 /iHz length, which were vertically 
stacked. The color bar indicates the power in cm 2 s -2 . The 
existence of two distinct groups of peaks separated by half 
the folding frequency is compatible with oscillation modes 
of consecutive degrees (e. g., 1 = 1 and 1 = 2). 
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Fig. 5. Top and bottom panels: zoom of the power spec- 
trum between 700 and 3400 /iHz. The dashed line indicates 
the background level fit and the vertical dash dotted lines 
indicate the frequencies of the 22 peaks with 1% significance 
level. 



compatible with the estimated thermal constant time of 
the Mach-Zehnder interferometer of approximately 30 min- 
utes. It appears that the 6000-/iHz region is modulated by 
a low frequency of 136 ± 5 /iHz, which corresponds to a 
period of about two hours. Such a periodic modulation of 
the guiding, equal to a fraction of the day, is commonly en- 
countered with telescope guiding systems and may explain 
the previo us incorrect determinatio n of the mean spacing 
frequency (jMosser et al.| [l993. 200QD . A similar modulation 
is not observed in the presumed p-mode region. In con- 
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Fig. 6. Autocorrelation of the power spectrum in three fre- 
quency ranges. The blue curve shows the low frequencies 
[10, 500] /iHz, which are dominated by Jovian rotation, 
and optical path difference fluctuations of the interferome- 
ter, generated by temperature fluctuations. The black curve 
shows the frequency range dominated by the pointing er- 
rors [4450, 7350] /iHz. The red curve shows the [900, 3800] 
//Hz range, where we measure a double excess power. Only 
this last curve displays a modulation at about 160 ±5 /iHz, 
indicated with the grey region, as well as at its half value. 



trast, we clearly identify two maxima at 80 ± 5 /iHz and 
160 ± 5 /iHz, which represents an alternative identification 
of the large spacing A^o- This result excludes an instru- 
mental artifact generating the comb structure. 

3.2. The contribution of solar oscillations 

SYMPA measures the Doppler shift of solar spectral lines 
after their reflection by Jupiter's cloud decks. Solar oscilla- 
tions may therefore be detectable in the data, particularly 
in the secondary excess power, given that solar p- modes 
extend from 1650 /iHz to 5000 /iHz, with a f requency of 
maximum amplitude at about 3050 /iHz (e.g. iGellv et al.l 
2002). At any given point on Jupiter, the total Doppler 
shift is the sum of the radial velocities in the solar atmo- 
sphere integrated over the entire solar surface, of Jupiter 
itself locally, and of the relative motions between Jupiter, 
the Sun, and the observer. However, because we consider 
the time series associated with antisymmetric modes, all of 
the global velocity fields (i.e. Jupiter's rotation and signa- 
tures from solar oscillations) should cancel out. 

To test the possibility of some leakage of the solar 
modes, it would have been useful to compare the so- 
lar and the SYMPA spectra, using the solar full disk 
Doppler measurements from the GOLF/SOHO instrument. 
Unfortunately, a four-day interruption in the GOLF data 
during the SYMPA run prevents us from directly compar- 
ing the two time series. BiSON data also suffered a 1-day 
interruption on April 3 2005 EL We however computed the 
spectrum of the GOLF data, for a period of 9.27 days pre- 
ceding the data interruption, with the same window func- 
tion. Although some of the peaks match solar frequencies 
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Fig. 7. Power spectra of SYMPA (top) and GOLF/SOHO 
data (bottom). The Jupiter data is a zoom of the power 
spectrum plotted in Fig. [2j The solar spectrum was ob- 
tained by using the GOLF/SOHO data on a 9.27-day run, 
in early April 2005, with the same window function and the 
same noise level as SYMPA: the GOLF data were interpo- 
lated on the same time grid, then a white noise with stan- 
dard deviation corresponding to those of SYMPA's data 
(18.6 m s _1 ) was added. 



to within ±10 /iHz, this behavior could be serendipitous 
and does not prove the presence of solar p-modes in the 
Jovian spectrum. Solar p-modes would be expected to be 
visible in the symmetric modes, and in particular in the 
time series corresponding to 1=0. However, this time series 
is affected by guiding and temperature fluctuations and ex- 
hibits a higher noise level in the [1000, 3000 /iHz] frequency 
range. Solar low-degree p-modes have maximum amplitudes 
of 20 cm s _1 and would be seen with half the sensitivity 
that we have for Jovian modes. Hence, owing our ~40 cm 
s _1 noise level for this range of frequencies, they are un- 
likely to explain the excess power measured between 2500 
and 3400 /iHz. In any case, they should not be considered 
responsible for the main excess power between 800 and 2100 
/iHz. 

3.3. The contribution from Jupiter's meteorological activity 

Since the detected signal is neither an instrumental nor a 
solar artifact, it appears to be of Jovian origin, as proven by 
the presence of twice the Jovian rotation in the comb-like 
structure (Fig. 2). However, Jupiter is a complex and dy- 
namic planet, with many sources of potentially misleading 
signals. The inhomogeneities in both the photometric and 
dynamic structure could indeed affect the measured signal. 
The effect of quasi-periodic albedo features and its possi- 
ble contamination o f Doppler measurements was studied in 
iLederer et al.l (|l995l ). However, they considered the effect of 
photometric pseudo-periodic features as possible source of 
spurious velocity signals in the case of spatially unres olved 
Doppler measurem ents as in ISchmider et al.l (|l99ll ) and 
M osser et al.l (|1993f ) . Most of the detected frequencies were 
found below 1000 /iHz, but that might be a consequence of 
the limited spatial resolution of their dataset. Unresolved 
features on Jupiter may indeed affect the Doppler signal 
and could possibly account for an appreciable amount of 
the high frequency signal currently identified as p-modes. 
These effects should not affect the signal we detect because 
the velocity map is computed before summation over all 



pixels and photometric contribution is thus cancelled out. 
As we simultaneously recorded the photometric images, it 
is possible to compute the power spectrum of the photo- 
metric time series associated with our measurements. For 
each exposure, the sum of the four outputs of the instru- 
ment is an image of Jupiter w ith the fringe pattern sup- 
pressed ([Schmider et al.ll2007l ). We multiplied each image 
for the whole run with the pair of masks l,m = 1,0 and 
l,m = 1,1, eliminating the outer 25% of Jupiter's diame- 
ter, to obtain two new times series. The power spectrum 
of the mean time series was then calculated with the dis- 
crete Fourier transform. Neither the power spectrum nor 
the echelle diagram associated with A^o present any simi- 
larity to the power spectrum of the radial velocity time se- 
ries (Fig. [8]). As expected, the peaks due to albedo features 
are integer multiples (48, 78, 82, for the most important) 
of the Jovian rotation frequency, corresponding to spatial 
pseudo-periodic features on Jupiter. The main peak in the 
photometric power spectrum (~ 2250 /iHz) corresponds to 
a range of minimum energy in the velocity spectrum. This 
proves that our processing fully suppresses the photometric 
fluctuations in the velocity measurement. Moreover, in the 
photometric signal, the absence of the comb-like structure 
detected in the velocity spectrum proves that the veloc- 
ity signal cannot be associated with the Jovian dynami- 
cal inhomogeneities that are associated with the pseudo- 
periodic albedo features. This is unsurprising as these dy- 
namical patterns have mainly horizontal components and 
would contribute to the Doppler measurement only at the 
limb, which was removed during the data processing. We 
therefore conclude that acoustic eigenmodes are the most 
probable explanation of the measured signal. 

3.4. Comparison with Jupiter's model 

The inferred Az^o value agrees with the predictions of mod- 
ern m odels for the interior of Jupiter ([Saumon fc Guillotl 
2004). This result is independent of the equation of state of 
hydrogen used for the modeling and even extends to mod- 
els that do not fit Jupiter's J4 measurement (Fig. [9]). Older 
interior models that used an equivale nt helium fraction to 
mimic the presence of heavy elements ([Gudkova et al.lll995l ) 
yield values of Auq that are slightly too small, thus showing 
that the present determination is already useful for validat- 
ing Jupiter models. The theoretical Az^o value depends only 
weakly on factors such as the core mass, but we note that 
the core composition parameter yields a rather significant 
~ 1/iHz change when varied between pure rock and pure 
ice models, which slightly favors the former solutions. More 
accurate determinations of A^o can thus provide with in- 
formation about Jupiter's central core, a crucial element to 
understand the planet's formation. However, the tightest 
constraints will be provided by the comparison of observed 
with theoretical individual mode frequencies. 

4. Discussion and conclusion 

4.1. Detection of Jovian p-modes between 800 and 2100 fiHz 

We applied the most sophisticated analysis methods of as- 
teroseismology to reanalyze the Jovian data, obtained by 
the SYMPA instrument, in April 2005, at Teide observa- 
tory. This led to the detection of a clear excess power be- 
tween 800 and 2100 /iHz and a secondary excess power be- 
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Fig. 8. Top panel: power spectrum of the photometric time 
series described in 2.2. As in Fig. [2j the thin grey and bold 
black lines represent the power spectrum smoothed with a 
3.1-/iHz wide triangular window, and a linearly increasing 
width triangular window, from 14 to 129 /iHz across the 
[500, 4500] /iHz frequency range, respectively. The spec- 
trum is normalized such that its maximum value equals 
1. The red dashed line indicates the background level fit. 
Bottom panel: echelle diagram of the power spectrum for a 
folding frequency of Az^o = 155 /iHz, calculated in the same 
way as in Fig. [5] 



tween 2400 and 3400 /iHz, which are both modulated by a 
comb-like structure. We could easily rule out the temper- 
ature and pointing effects, and albedo features on Jupiter 
as causing the measured features in the power spectrum 
of the photometric time series associated with the dataset. 
Solar-p modes are similarly unable to account for the main 
excess power. 

In terms of the secondary excess power, the short ob- 
serving time associated with a low duty cycle, prevents us 
from being able to explain such a double bump in terms of 
the three following hypotheses: 1) It is a snapshot of the 
modes that were stochastically excited in early April 2005. 
2) It reveals a variation in the mode amplitude as a func- 
tion of frequency. This hypothesis cannot be excluded in the 
case of Jupiter, as the inner structure could contain several 
discontinuities that could affect the quality factor of modes 
with different radial order. 3) The secondary excess power 
is the reflection of solar modes. The last point is the least 
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Fig. 9. Masses of the central core and values of A^o pre- 
dicted by interior models of Jupiter. Colored symbols repre- 
sent solutions obtained for different EOSs of hydrogen (as 
labeled) and using consisten t EOS s for helium, ices, and 
rocks (see lSaumon fc Guillotl (|2004l ) for a full description of 
the models). Boxes, open squares, and open circles corre- 
spond to solutions matching all available constraints. Dots 
correspond to models that fail to match the observed fourth 
order gravitational moment J4. The models assume a two- 
layer hydrogen-helium envelope overlying a core made ei- 
ther of pure rocks (circles) or pure ices (squares). Asterisks 
correspond to older model solutions i n which ices and rock s 
were modeled using the helium EOS ([Gudkova et al.l lT995). 
The grey area indicates the 1-a constraint obtained from 
our determination of the large separation of the modes. 



probable because the solar large-spacing signature does not 
appear in our data. Moreover, since the cut-off frequency 
of Jo vian modes is about 3500 /iHz ([Lognonne fc Johnson! 
120071 ), finding modes at these frequencies would not be sur- 
prising. Finally, the secondary excess power also mirrors the 
double ridge structure that is observed in the main excess 
power in the echelle diagram (Fig. [5j). 

In any case, the most probable hypothesis is that 
the signal between 800 and 2100 /iHz is generated by 
Jupiter's acoustic modes. We measured the global seis- 
mic parameters, i. e., the frequency of maximum ampli- 
tude at ^max = 1213 ± 50 /iHz and the mean large spacing 
A^o = 155.3 ±2.2 /iHz. The large spacing value agrees 
with theoret ical models o f the internal structure of Jupiter. 
In addition, Mosser (1995) calculated the quality factor of 
the Jupiter cavity, then infered that the mode lifetime is 
longer than 12 days. Thus, the mode damping is almost 
negligible during a nine-day run and we estimate the mode 
amplitude using the pure sine assumption, which causes the 
maximum mode amplitude to be 49^ cm s_1 - The SNR 
of the detection is similar to the first asteroseismic results 
for solar-like targets, from data acquired by the CoRoT 
satellite (iBaglin et al.ll2002[) ; in pa rticular those obtained 
for HD 181906 (jGarcfa et al.ll2009l ). 

The large error bar in the peak frequencies (—5/ + 
9 /iHz) prevents us from identifying the modes. From Table 
[TJ the leakage between the masks implies that the de- 
tected peaks could correspond to modes of degree (1 = 1, 
m =±1), (1 2, m=-l,-2), (1 = 3, m=0, 1,2,3), and (1 
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= 4, m=— 4, —3,-2,-1), which are almost uniformly dis- 
tributed as a function of frequency, by considering that the 
Jovian rotation introduces a frequency splitting of 28 /iHz 
and that modes of consecutive degrees are spaced by about 
half the large spacing. 

4.2. A new start for giant planet seismology 

To confirm the present observations, further observations 
are planned with an improved model of SYMPA, including 
an accurate pointing system, for the 2012 and 2013 Jovian 
oppositions. However, seismic measurements will strongly 
constrain Jupiter's inner density profile only with data of a 
significantly higher SNR, when modes are unambiguously 
identified in the spherical harmonics base. The main limi- 
tation of ground-based measurement arises from the poor 
duty cycle and atmospheric turbulence. Obviously, the op- 
timal observing conditions would be obtained from space 
observations, which is th e purpose of the Doppler spectro- 
imag er DSI/ECHOES (ISchmider fc DSI /ECHOES Team 



2010 ) proposed to the ESA Laplace mission ([Blanc et al 
2009h . 



In any case, our results provide the clearest evidence 
yet of Jovian global oscillations. These measurements open 
the possibility of directly probing the planet's interior and 
determine its global composition, core size, presence of 
phase transitions, and even internal rotation. The tech- 
nique should be applicable to all giant planets in the Solar 
System, with unique prospects for the understanding of 
their composition and formation. 
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